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I. INTRODUCTION 

In our four-dimensional universe 'real' astrophysical 
black holes cannot carry large electric charge. In the 
presence of surrounding plasma, flux of the plasma par- 
ticles with the opposite charge rapidly reduces the orig- 
inal charge of the black hole and makes it small. In the 
absence of monopoles a black hole cannot also carry a 
magnetic charge. However, in the presence of accret- 
ing disk, magnetic held can exist in the black hole vicin- 
ity. Moreover it is plausible that it plays an important 
role in a black hole activity. The electric and magnetic 
fields E and B carry energy density e ~ \E\ 2 + \B\ 2 , 
but, in a realistic situation its contribution to the curva- 
ture Ge/c 4 is much smaller than the characteristic curva- 
ture of the spacetime near the horizon ^(grav. radius) -2 . 
This means that the distortion of the background gravi- 
tational field is small and can be neglected. We call such 
electric and magnetic fields weak. This does not mean at 
all that their role in the physical processes in the black 
hole vicinity is also negligibly small. The reason is that 
for charged particles the relative weakness of the electro- 
magnetic field can be compensated by a large value of the 
chargc-to-mass ratio, which, for example, for electrons 
is e/m w 5.2728 x 10 17 g -1 / 2 cm 3 / 2 s -1 . In other words, 



the motion of charged particles in weakly charged and/or 
weakly magnetized black holes might be quite different 
from the motion of neutral ones. 

For example, the innermost stable orbit of a charged 
particle moving in the equatorial plane of a weakly mag- 
netized black hole can be very close to the horizon, see, 
e.g., 0. In consequence, as it was shown recently, weakly 
magnetized black holes may play a role of particle accel- 
erators dm. 

The study of weakly charged and/or magnetized black 
hole is greatly simplified, since the electromagnetic field 
can be considered as a test one. The electromagnetic field 
A in the Lorentz gauge V a A a = obeys the equation 



UA a - Ric 6 Q A b = . 



(1.1) 



If the spacetime possesses a Killing vector £, it is easy to 
show that it satisfies the equation 



□f„ + Ric 6 Q & = . 



(1.2) 
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Hence in a Ricci-flat spacetime £ is a potential for a spe- 
cial test electromagnetic field. This observation by [aQ 
is well known and often used. Let us emphasize that this 
property is valid in a spacetime with arbitrary number of 
dimensions. A well known example is the Kerr spacetime 
with two Killing vectors, £( t ) and generating time 
translation and axial rotation symmetries. The corre- 
sponding test Maxwell fields describe respectively a weak 
electric field of a test charge inside the black hole, and 
the axisymmetric stationary magnetic field near the black 
hole which is homogeneous at the infinity and directed 
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along the axis of rotation. Evidently, this construction 
can be easily generalized to the Myers-Perry spacetimes 
in an arbitrary number of dimensions, e.g., [7Hl0j. 

In spite of the fact that the relation between special 
test fields and the Killing vectors generating them is of a 
general nature, the character of the motion of a charged 
particle in the presence of such a field is quite differ- 
ent for different fields. To illustrate this, let us consider 
again the case of 4D Kerr spacetime. The geodesic equa- 
tions for particles in the Kerr metric are completely inte- 
grable llj. This property is preserved for the motion of 
electrically charged particles in a weakly charged black 
hole. (In fact, it is valid for the Kerr-Newman space- 
time, where the charge of a black hole is not weak.) The 
case of a weakly magnetized Kerr black hole is quite dif- 
ferent. Particles moving in the equatorial plane are still 
completely integrable, while the 'fourth' (Carter's) inte- 
gral of motion is absent for the generic motion out of 
the symmetry plane. As a result, in a general case the 
motion of a charged particle in weakly magnetized black 
holes is chaotic. 

It is interesting, that this is a generic feature of higher 
dimensional rotating black holes as well. Namely, it was 
shown in [l[ that the charged particle equations in higher 
dimensional rotating black holes with a test electric field 
proportional to the primary (timelike) Killing vector field 
are completely integrable — the properties of such space- 
times are similar to the properties of a weakly charged 
Kerr spacetime. However, complete integrability is de- 
stroyed in the presence of a test electromagnetic field 
generated by secondary (rotational) Killing vector fields. 
(This is an analogue of weakly magnetized black holes.) 
It was shown that the Klein-Gordon and Hamilton- 
Jacobi equations for such a system allow complete sepa- 
ration of variables only for the primary Killing test field 
[l[ . The latter stems from the existence of the Lax tensor 
for a charged particle — discussed recently in fl2| . 

In this paper we study the properties of a charged 
Dirac equation in higher dimensional rotating black hole 
spacetimes in the presence of a test electromagnetic field 
derived from various Killing vector fields, and demon- 
strate, that the structure outlined above for particles and 
scalars remains similar. Namely, for the weakly charged 
black holes the Dirac equation preserves its remarkable 
property and allows a complete separation of variables. 
Whereas we show that this is not possible for a magne- 
tized black hole. 

The paper is organized as follows. In the next sec- 
tion we review the charged Dirac equation in curved 
space and discuss its first order symmetry operators. In 
Sec. IIIII wc gather the information about weakly charged 
and magnetized rotating black holes in any dimension. In 
Sec. lIVI we demonstrate the existence of a complete set of 
first order mutually commuting operators for the charged 
Dirac equation in a weakly charged black hole spacetime. 
The explicit representation of these operators is found 
in Sec. [V] Separability of the charged Dirac equation in 
weakly charged black hole spacetimes is demonstrated in 



Sec. IVI1 We conclude in Sec. VII. Appendix [Al displays a 
derivation of the important identity which is crucial for 
showing the mutual commutation of the symmetry oper- 
ators in the set. Appendix [Bl offers an alternative proof 
of this property employing the explicit representation of 
these operators. In many places the paper closely follows 
fl3j where similar results were proved for the standard 
(uncharged) Dirac equation in higher-dimensional Kerr- 
NUT-AdS spacetimes. 

II. CHARGED DIRAC EQUATION IN CURVED 
SPACE 

A. Dirac bundle 

In what follows we parametrize a dimension of the 
spacetime as 

n = 2W + £, (2.1) 

with e — 0, 1 corresponding to the even, odd dimension, 
respectively. The Dirac bundle DM has fiber dimension 
2 N . It is connected with the tangent bundle TM of the 
spacetime manifold M through the abstract gamma ma- 
trices 7 Q € TM ®D\M, which satisfy 

7 a 7 b + 7 b 7 a = 2 g ab ■ (2.2) 

They generate an irreducible representation of the ab- 
stract Clifford algebra on the Dirac bundle. All lin- 
ear combinations of products of the abstract gamma 
matrices (with spacetime indices contracted) form the 
Clifford bundle CIM, which is thus identified with the 
space D\M of all linear operators on the Dirac bundle. 
The Clifford multiplication ('matrix multiplication') is 
denoted by juxtaposition of the Clifford objects. The 
gamma matrices also provide the Clifford map 7* , which 
is an isomorphism between the exterior algebra AM and 
the Clifford bundle CIM, 

if, = 7 *w = ]T i K) ai ... ap7 Q1 - Qp • (2.3) 
p 

Here, lo = ^ p uj p G AM is an inhomogeneous form, uj p its 
p-form part, lo p € A P M, and 7 ai --- a p = -y^ 1 ■ ■ ■ 7°?]. For 
future use we also define an operator 7r as -klo = ^2 p P^ P - 
We denote the standard (uncharged) Dirac operator 
both in the exterior bundle and the Dirac bundle as D 

D = e a V a , D = 7,e Q V a = 7 a V a . (2.4) 

Here e a € TM ® AM is a counterpart of 7 a in the exte- 
rior algebra and V denotes the spinorial covariant deriva- 
tive. We shall also denote by X a the object dual to e a 
and by A a contracted wedge product, defined inductively 

as [HQl 

aA/3 = aA/3, a A /3 = (A a _i a) A (A Q _i /3) , (2.5) 

k k — 1 
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where the 'hook' operator j corresponds to the inner 
derivative. The reader is also referred to the Appendix 
of UM for more details on the notation. 



where <p is an arbitrary function satisfying the 1 anomaly 
condition' 



dip + {dA) A K = , 



(2.11) 



B. Killing— Yano tensors 

A Killing-Yano (KY) tensor / [l?} is a p-form on the 
spacetime, the covariant derivative of which is completely 
determined by its antisymmetric part, i.e., by its exterior 
derivative df: 



^afai...a r ^[a/a 1 ...o r ] p+1 



Jaai ...a p 



(2.6) 



A closed conformal Killing-Yano (CCKY) tensor h is 
a p-form on the spacetime, the covariant derivative of 
which is entirely determined by its divergence £: 



p — I V h a 

<5a 2 ...a r n — r+l Va a 2 ...a r 



(2.7) 



KY and CCKY tensors are related to each other through 
the Hodge duality: the Hodge dual of a KY form is a 
CCKY form and vice versa. As we shall see, both these 
tensors play a crucial role in the study of symmetry op- 
erators of the Dirac operator. 

It is also useful to introduce the notion of a principal 
Killing-Yano (PKY) tensor. We define this to be a non- 
degenerate (with maximal matrix rank and functionally 
independent eigenvalues) rank-2 CCKY tensor [l8j]. We 
shall see in Sec. IIIII that this tensor exists in rotating 
black hole spacetimes in all dimensions. It generates a 
Killing vector £, 



1 



n- 1 



(2. 



which we call the primary Killing vector, as well as, a 
whole set of secondary Killing vectors, see Sec. IIIII 



C. Symmetry operators 

We shall now write down the most general first-order 
operators commuting with the U(l) charged Dirac opera- 
tor. The results follow from the general discussion in [l5[ 
where the symmetry operators of a Dirac operator cou- 
pled to an arbitrary flux were studied (for special cases 
see also [lj| and (20|). 

Let V be the Dirac operator for a charged particle in 
the U(l) potential A, 



V = D-iA = e a V a , V = j a V a 



(2.9) 



where we have introduced the U(l) covariant spinorial 
derivative V a = V a — iA a . 

We first consider the case when the background met- 
ric admits a Killing vector equivalent to a Killing-Yano 
1-form k. Then, the following operator: 

K K = X a -\nV a + jdK - iAAK + iip, (2.10) 



commutes with the Dirac operator T>, (|2.9[) . Given the 
potential A, one may try to find such ip so that the 
anomaly vanishes. Two possible choices are obvious 
(though not necessary): i) One chooses cp = in which 
case the operator and the anomaly condition become 



K K =X a AnV a ^ 
= k A dA . 



1 



-dn . 



(2.12) 



We shall see that a similar condition is required for 
higher-rank Killing-Yano tensors as well, ii) One chooses 
ip = AAk. In this case the operator (|2.10p and the 

anomaly (|2.1ip simplify to 



K K = A a _i nV a + ^d K . 
= £ Kl A, 



(2.13) 



where £ K « is the Lie derivative along the Killing vector 
k*. Hence, if A respects the symmetry of k, the opera- 
tor (|2.13[) defines a symmetry operator (i.e., an operator 
commuting with T>) for the charged Dirac equation. In 
particular, as we shall see in Sec. IIIII this property is sat- 
isfied in weakly charged and magnetized rotating black 
hole spacetimes for A associated with any Killing field, 
as all of them do mutually commute. 

Since any Killing vector is Hodge dual to an 
(n— l)-rank CCKY tensor we do not have to discuss the 
case of CCKY tensors of this rank anymore. 

Having dealt with these special cases separately, we 
may now formulate the following result: The most gen- 
eral (different from that discussed above) first-order op- 
erator S which commutes with the Dirac operator D, 
[D, S] = 0, splits into the following Clifford even and odd 
parts: 



where 



S c = L fo =X a _if V a 



7T - 1 



dfo, 



(2.14) 



(2.15) 



2tt 

S = M he = e a A ft, V a - ? '~ 7r ~, 1 5h c , (2.16) 

2{n — 7r) 

with / being an inhomogeneous odd (at least rank 3) 
KY form, and h e being an inhomogeneous even (at most 
rank- (n— 2)) CCKY form, subject to the condition 

(cL4)A/ o = 0, (dA)Ah e = 0. (2.17) 



Note that in odd dimensions the operators L and M are 
related by the Hodge duality. However, this is not true 
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in even dimensions, where these operators are necessarily 
independent. 

Indeed, in an odd number of spacetime dimensions 
the Hodge duality of Killing- Yano tensors translates into 
the corresponding relation of symmetry operators L and 
M. Namely, let z be the Levi-Civita rt-form satisfy- 
ing z ai ,,, an z ai " ,an = n!0 Then, in odd dimensions, the 
Hodge dual of a p-form w can be written as 

* u = zlj, (2.18) 

and for the operators of type L and M it holds that 



L zh = z M h , M x f = zLf, 



(2.19) 



where / is an odd KY form and h an even CCKY form. 
In particular, this is also valid for Killing vectors and 

dual CCKY (n-l)-forms. 

The anomaly conditions (|2.17|) are also related by the 
Hodge duality, thanks to the identity 

F A(*a) = *(F Aa) , (2.20) 
which holds for any 2-form F and an arbitrary form a. 



III. WEAKLY CHARGED ROTATING BLACK 
HOLE SPACETIMES 

In what follows we shall study the Dirac equation in 
the vicinity of the weakly charged and weakly magnetized 
rotating black hole spacetimes in all dimensions. Namely, 
as a background metric we shall use a 'canonical space- 
time ', which is the most general spacetime admitting the 
PKY tensor [2l|, [22| . When vacuum Einstein equations 
are imposed such a spacetime describes the most general 
known rotating Kerr-NUT black hole discovered in [23l j. 
The metric can be equipped with a weak electromagnetic 
field associated with (N + e) isometries present in the 
spacetime. Among them the one associated with the pri- 
mary Killing vector plays a special role. The background 
metric with vector potential proportional to this Killing 
field describes a weakly charged rotating black hole. We 
shall see that in this case the charged Dirac equation 
separates and the separability can be characterized by 
the existence of a complete set of mutually commuting 
first order symmetry operators. On the other hand, for 
weakly magnetized black holes, which are described by 
a vector potential associated with the secondary Killing 
vectors, the first order operators do not mutually comm- 
inute and the Dirac equation cannot be (at least in the 
canonical frame) separated. It is very likely that in this 
case the integrability is lost and the Dirac equation does 
not separate in any frame. 



1 Note that 7* (z) is the ordered product of all n gamma matrices 
and in odd dimensions it is proportional to the unit matrix. See 
also formula (15-101 ). 



A. Canonical spacetimes 

The most general metric admitting the PKY tensor, 
also denoted as a canonical metric, writes in canonical 
coordinates as |2lL Y£M 



N 

where 



JV-l 



Ql 



3=0 



(3.1) 



3=0 



Ql 



U u 



4 fc) = E 



ui<—<u k , Vi^n 



, AW ^ 

v lt ,,,,v k 



s 



AW 



(3.2) 



Coordinates x^ (fi = 1, . . . , N) stand for the (Wick ro- 
tated) radial coordinate and longitudinal angles. Killing 
coordinates tpk (k = 0, . . . ,N — 1 + e) denote time and 
azimuthal angles associated with Killing vectors £/m 



(3.3) 



At points with x^ = x v with fj, ^ v the coordinates are 
degenerate. We assume a domain where x^ 7^ x v for 
/i^i/. In such a domain we can always order and rescale 
the coordinates in such a way that 



+ %v > and x u 



x u >0 for fx < v . (3.4) 



With this convention and assuming Euclidian signature 
we have 

^ = (-1)^1^1, X, = (-1) N -»\X,\. (3.5) 
The metric can be written in the diagonal form 
g = ( E» <8> E 1 * + £ A ® E» \ + e E° ® E a , (3.6) 

where we have introduced an orthonormal covector frame 
E a = {E^,E^,E }, 



N-l 



E° = y/sJ2^ {j) dipj ■ 
3 

The dual vector frame E a = {E^, Ep,, Eo} is 



(3.7) 



En 



3 



V,- 



(3.8) 
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Note that E° and Eg are defined only in an odd dimen- 
sion. The Ricci tensor is also diagonal in this frame and 
is given by [H| 

Ric = J2 r »( E>i ® E " + E(l ® E ^) + £r o E " ® E° , 



f 

2x~, 



E 


4( 


x v * 






x P ( 


E 





x £C 

II 



(3.9) 



The PKY tensor of the canonical metric reads 



h = E x>1 Efi A ^ A 



(3.10) 



It generates a tower of even rank CCKY forms and 
even or odd rank (in even or odd number of spacetime 
dimensions) KY tensors fV>, 



■ h A ■ • • A h 



(3.11) 
J- 

fU) = t h U) = (-1)3 zh^ , (3.12) 
as well as the 2nd-rank Killing tensors (j = 0, . . . , N — 1) 

kU) = E A » ] s " ® E ^+ E(l ® eP ) + £A^E° $ E° . 

_ (3.13) 

The primary Killing vector £, (|2.8|) . coincides with £( ) 
given by (|3.3|) . £ = £(o) = d^ , and can be explicitly 
written as 



i = E y/o* E n + s^sbo ■ 



It satisfies an important relation 
1 



n-2j+l 



(3.14) 



(3.15) 



The secondary Killing vectors can also be generated from 
the PKY tensor, e.g., [IH, and are explicitly given by 



(3.16) 



We conclude with a few identities that will be used in 
Sec. IIVI The following two relations have been proved in 

[hW,h®]=0, (3.19) 
[e (a A/i (fc) , e b) Ah (l) ] = 0. (3.20) 

The latter identity implies 

[lj A h^ k \ uj A h^] = (3.21) 

for any 1-form uj. In Appendix [A] we also prove 

V [o &] [e[ Q A/iW, e b lA/iW] + =0. (3.22) 

In all these identities, [ , ] and [ , ] + denote Clifford com- 
mutator and anti-commutator, respectively. Other prop- 
erties of canonical metrics are gathered, for example, in 
recent reviews 0, [26[ . 

B. Test electromagnetic field 

We shall now equip the canonical spacetime with the 
test electromagnetic field A given by a linear combination 
of background isometries £}- k \ 



A 



(3.23) 



with constants qk characterizing the field strength. In 
order this field to be a solution of the vacuum Maxwell 
equations, we have to impose a condition Ric a ;,yl'' = 0, 
cf. Eqs. (|1.1[) and (jl.2[) . Taking into account the diago- 
nal form (|3.9[) of the Ricci tensor, the explicit form (|3.16[) 
of the Killing vectors and a functional independency of 
the functions A^} , we find that spacetime must be a vac- 
uum Kerr-NUT spacetime [23| . Ric = 0, with the metric 
functions Xn in the forrri 



JV-l 



X, 



- 26p x^ 



l-s 



ec 



(3.24) 



k—e 



It describes a general rotating black hole with spherical 
horizon topology. Parameters Cfc , c and are related to 
(N+e) independent angular momenta, a mass parameter, 
and NUT charges. 



In odd dimensions we have 



cf 



(N) 



All the Killing vectors mutually Lie-bracket-commute, as 
well as preserve all the Killing- Yano and Killing tensors 
constructed from h, e.g., 



2 We impose the condition of vanishing electric current to guar- 
(3.17) antee a realistic matter content of the spacetime. However, the 

commutativity of the symmetry operators and the separability 
of the Dirac equation discussed below are independent of this 
requirement. For example, we can relax the vacuum condition 
for the metric and consider a nontrivial cosmological constant. 
In such a case, the spacetime must be filled with a charged test 
fluid, with electric current aligned with the vector potential, i.e., 
(3 18) with the chosen isometry of the spacetime. 
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After fixing the metric, let us now study what restric- 
tions the test electromagnetic field A has to satisfy in 
oder it to be compatible with the existence of n first or- 
der operators commuting with the charged Dirac opera- 
tor. Applying the results of Sec. [TTl in any dimension the 
candidates for such operators are generated by (N + e) 
Killing forms £(*) and N CCKY tensors E The corre- 
sponding operators K^ k ) , cf. (|2.13j) . and M h{j) , cf. (|2.16l) . 
commute with T> only if the anomalous conditions (|2 . 13[) 
and (|2.17p are satisfied. 

Since all the isometries ^ha of the Kerr-NUT spacetime 
mutually commute [first property (|3.18[) ]. the anomaly 
condition 

C kk) A = Q (3.25) 

is satisfied and the operators (|2.13[) commute with T> for 
any A given by (|3.23[) . 

However, for the operators (|2.16|) the situation is more 
complicated. Let us start from the second property (|3.18[) 
and re-express the Lie derivative of h along A as follows: 

= C A »h = d(AAh) 

= —dAA h + V a A A A a _i h + V A »h 

i 

= -^dAAh + AA^, (3.26) 

where we have used the fact that A is a Killing form 
and h a PKY tensor. Obviously, only for A proportional 
to the primary Killing form £ b , the latter term vanishes 
and the anomaly condition dAA h = is satisfied. The 

same remains true for higher rank CCKY tensors h^'. 
Namely, the anomaly condition 

dA/\h( j) =0 (3.27) 

is satisfied only for such A given by (|3.23|) which is pro- 
portional to the primary Killing form, 

A = e£, b , (3.28) 

i.e., Qo = e and Qk = otherwise. Hence, from the 'math- 
ematical point of view' the test electromagnetic field 
determined from the primary Killing vector is funda- 
mentally different from the test electromagnetic field 
determined by the secondary Killing fields; only when 
A = e £ , all n first order operators (|2.13[) and (|2.16l) com- 
mute with the charged Dirac operator T>. 

Physically, the test electromagnetic field (|3.28p given 
by the primary Killing vector describes a weakly charged 
Kerr-NUT spacetime, whereas the test field associated 
with the secondary Killing vectors corresponds to a black 
hole immersed in a magnetic field. 



3 Odd KY tensors f^> in odd dimensions are dual to even CCKY 
tensors few. Hence due to the fact that both, the operators 
112,191 1 and the anomaly condition I I2.17I I. are Hodge dual, with- 
out loss of generality we may focus on operators generated from 
CCKY tensors only. 



C. Charged scalar particle 

It is well known that in the canonical spacetime in the 
absence of a test electromagnetic field, the Hamilton- 
Jacobi equation as well as the Klein-Gordon and Dirac 
field equations allow separation of variables [13, [28| ■ 

The motion of a scalar charged particle in the weakly 
charged or magnetized spacetimes has been studied more 
recently It was observed that in the weakly charged 
Kerr-NUT spacetime the motion of a charged scalar par- 
ticle maintains its complete intcgrability, while the cor- 
responding charged Hamilton- Jacobi equation separates. 
This integrability can be characterized by the existence of 
a nontrivial Lax tensor discussed in [12| . The separabil- 
ity property remains true also for a charged scalar field 
in the vicinity of the weakly charged black hole. Such 
separability is underlain by the existence of a complete 
set of mutually commuting symmetry operators gener- 
ated from Killing vectors and Killing tensors 
§T3]), 

= - [Va - iA a ) k^ ab [V 6 - iA b ) , (3.29) 

with A given by ([3~28]l . 

On the other hand, it has been observed in [l| that sim- 
ilar operators do not commute when instead of the pri- 
mary Killing field ^ , A is associated with the secondary 
Killing fields, i.e., in the immersed magnetic field. In 
this case neither the charged Hamilton- Jacobi, nor the 
charged scalar field equations separate in canonical coor- 
dinates. 

The aim of this paper is to demonstrate that a similar 
conclusion remains valid for a charged spin ^ field. We 
already know that in the weakly charged case there exist 
ri—1 operators commuting with the charged Dirac oper- 
ator T>. In the next section we demonstrate that all these 
operators mutually commute and hence together with T> 
form a complete set of symmetry operators. In Sec. IVII 
the explicit separation of the charged Dirac equation in 
the chosen representation will be demonstrated. 



IV. COMPLETE SET OF DIRAC SYMMETRY 
OPERATORS 

A. Complete set of operators 

The canonical spacetime (|3.ip in the absence of a test 
electromagnetic field allows a complete separation of vari- 
ables for the Dirac field equation [28|. It was shown in 
[HI, [l6| that such a separability can be characterized by 
the existence of a complete set of mutually commuting 
first order symmetry operators. These operators are de- 
termined by the tower of symmetries built from the PKY 
tensor h. Namely, they are given by (N + e) KY 1-forms 
C {k \ (O, and N CCKY forms ((3~TT]) . In the exte- 
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rior algebra notation they read: 

K^=K^=X^^V a + \d^, (4.1) 
Mj 0) = M<°>, = e Q A h«)Va + \{n-2j) £ A fc^-i) 

In the expression for Mj 0) we have employed (|3.15|) to 
simplify the second term. Note also that the operator 
A/q ^ coincides with the Dirac operator, Mo = D. 

In the case of weakly charged Kerr-NUT spacetimes 
discussed in the previous section, the operators (|4.ip can 
be upgraded, to give a complete set of mutually com- 
muting operators of the charged Dirac operator V, (|2.9I) . 
Following the results (gH5J| and (|2?ToTl from Sec. Hfl to- 
gether with the discussion in Sec. IIII B[ the operators 
commuting with the charged Dirac operator are 



K k = A>, = Kf ] , 
Mj=M h(J ) =Mj 0) -ief AfcW , 



(4.2) 
(4.3) 



with k = 0, . . . , JV - 1 + £ and j = 0, . . . , JV — 1. We 
shall now show that, similarly to the uncharged case, the 
operators (|4.2|) and (|4.3|) mutually commute. 



B. Mutual commutativity 

It was shown in [la | that the commutator with the op- 
erator if K is equivalent to a Lie derivative along a Killing 
vector k". The mutual commutativity of (|4.2[) thus fol- 
lows from the Lic-commutativity of the Killing vectors, 
first property (|3.18[) . while the commutativity of opera- 
tors (|4.2|) with operators (|4.3[) follows from the second 
property (|3.18j) . 

It remains to prove that also all the operators (j4.3[) 
mutually commute. The commutator can be split into 
the following six terms: 

[M k , M t ] = [M< 0) , M/ 0) ] 

-2ze<£ a [e (a A/i (fc) , e b) A h^] V fc 
+ 2ze(V [Q 4])[e [a A/ i W, e b l A^] + 

n-2fe ri . , , t _n . . , mi (4.4) 
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-e 2 [^AfcW.fA^], 

where we have used £ A V a /i (fe) = 0, cf. (|3~T5|) . Employing 

the commutativity of uncharged operators , and the 
identities (|3.20p ~ p.22p . we find that all the terms inde- 
pendently vanish and hence the operators (|4.3p commute. 
[An alternative proof using the explicit representation of 
operators Mj is gathered in App. [5]] 



V. EXPLICIT REPRESENTATION 

In this section we find an explicit representation of the 
action of operators (|4.2p and f)4.3[) on the Dirac bundle 
in which the charged Dirac equation separates. As cus- 
tomary, we shall omit the explicit symbol 7*, defined by 
(|2.3p . whenever it is clear from context that operators 
are being considered instead of forms. So, for example, 
we write Kj for j*Kj and Mj for 7*Mj. 



A. Spinors and 7 matrices 

Following [l3| , we represent the fiber of the Dirac bun- 
dle as a tensor product of N 2-dimcnsional spaces S, 
DM = S N M. We use Greek letters e = ±, ? = ±, . . . for 
tensor indices in these 2-dimensional spaces and write a 
generic 2-dimensional spinor as x — X + $+ + X~$- = 
X e i9 € . Here and i?_ represent a frame in the 2- 
dimensional spinor space S and the components of the 

spinor are two complex numbers ( V The Dirac spinor 

ip e DM is then written as 



V = ^ ei ™ e *0«... ejv (5.1) 

with 2 N components ip ej ---- e N and the Dirac bundle frame 
is given by 



# ei (8 • • • ® i? e 



(5.2) 



The gamma matrices are constructed as various tensor 
products of Pauli matrices. Let 1, t, cr, and a be the unit 
and Pauli operators on S, respectively. In components, 
their action reads 



(IX) £ =X £ , (ix) £ =ex £ , 
{o xY = X~ £ , (o- xY = -ie X" 
In matrix form they are written as: 



(5.3) 



1 

1 

1 

1 



1 

-1 

— i 

1 



(5.4) 



For any linear operator a G S\M we denote by 
ot/fj,) S D\M a linear operator on the Dirac bundle 



1 ® a <£> 1 1 



(5.5) 



with a on the /_t-th place in the tensor product. Similarly, 
for mutually different indices . . . , jij we define 



(5.6) 



Using these preliminaries, the abstract gamma matri- 
ces with respect to the frame £7 a = {E^, Ej±,Eo} can be 
written as 



7 M = '-<i... m -i> (J (m> > ^ = L {i...n-i)°(n) , 
7° = 1{1...N) , 



(5.7) 
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with 7 defined only in odd dimension. In components, 
the action of these matrices on a spinor (|5.ip reads 



u-l 



( 7 "V0 61 - 



(7^) £1 - 



.(-e A )...ew 



/i-l 



.(-e fI )...e JV 



(5.8) 



N 



i/=i 

Note the relations 

and the fact that in odd dimensions we have 

7,(z) = yi-PArAi-Aw o = f JV _ ( 5 10 ) 



(5.9) 



B. Symmetry operators 

Symmetry operators (|4.2p corresponding to Killing 
vectors are, in general, equivalent to the Lie derivative 
lifted from the tangent bundle to the Clifford or Dirac 
bundles. Thanks to (13.31) we thus write 



K k 



a 

dipk 



(5.11) 



where is a partial derivative along ip k , acting just on 
the components of the spinor with respect to frame #e, 
i-e., K kX = JL- X = %L* B . 

To lift the operators Mj, (|4.3|) . to the Dirac bundle 
we employ (|2.3|) . Using the results in [l3| we have the 
following important formula valid for any 1-form a: 



(5.12) 

Here, B^ and fiW represent a 'spinorial analogue' of 
functions and ((372]) . and read 



E 



l '(v- i ) x vi ' ' ' L (v k ) x i>k i 



"1, 

Ul< — <U k , Vi^tl 



E 



l {v\) x v\ • • • L (u k ) x Vk 



(5.13) 



"1,— 

I/l<"'<ffe 



In particular, we have 
7* (V A/i ( ^ 



(5.14) 



,7+1 



1 ( Ev^^'^w/ + ieVSB^ 



Next, one has to use the expression for a spinorial deriva- 
tive 



V a = 5 a + |W abc 7 f, 7 c , 



(5.15) 



where 3 a is a derivative acting only on components of the 
spinor and the spin connection coefficients are gathered 
in appendix of [13]. The resulting formula for Mj reads 



x 'n L i 



i 



4X„ 2 



J £t ^^v) x v 



II 



+£ V 



2 3 Ett: 



(5.16) 



.7 + 1 



gO'-i) _ g(J) f 



ic dip 



-2e 



VI. SEPARABILITY OF THE CHARGED 
DIRAC EQUATION 

Let us formulate the main result: the commuting 
symmetry operators and Mj have common spinorial 
eigenfunctions tp 

K k ^ = i^ k ^, (6.1) 
MjV = Xjip , (6.2) 
which can be found in the tensorial R-separated form 

V = R exp(i J2 k ■ (6.3) 

Here is an iV-tuple of 2-dimcnsional spinors, where 
each spinor Xv depends only on the variable x v , Xu = 
Xu{x u ), and R is the (Clifford bundle)-valued prefactor 



R 



n( 

K.X 



(6.4) 



Eqs. (|6.ip and f|6.2[) are satisfied if and only if the spinors 
Xv satisfy the ordinary differential equations (|6.20l) be- 
low, which generalize the conditions found in [28j and 
fl3j for the case e = 0. 

The spinor ip given by (|6.3[) clearly satisfies (|6.1j) . To 
show (|6.2p . we need to calculate Mjtp, with Mj given by 
([576]) . We have 



MjV' 



1 ex p(* Y^k^kipi 



(6.5) 



L M 



x f 3 i 3e i 1 y 1 i *" 



iy/S 



i g(j- 1 )_/3(i) 



2* 



N 



7' 



RQQx» 



where we have performed the derivative with respect to 
angles %p k and introduced the functions ^ M of one variable 
x M given by 



^ = E^(-^)™'-^ 



(6.6) 
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Let us now concentrate on the derivatives of the prefactor 
R. Using the following relations proved in (l3j : 



r 



where 



dx. 



V„, = 



It. 



2 Xp, i^u)X u 



we can bring the operator R to the front, to get 



(6.7) 
(6.8) 

(6.9) 



E^-^f"^ 



( 



(6.10) 



Z ^ y 



e 7 



which is to be compared with 

Xjtp = Xj cxp(i J2 k ^kipk)R( 



(6.11) 



Following jl 31 ] we further introduce the following func- 
tions X v of a single variable x v : 



(6.12) 



Note that in odd dimensions the constant Xn, defined 
by Mj^tp = Xj^ip, is not independent. In fact, using the 
relation between operators K and L, the Hodge duality 
(|2TT9"|) together with (pTTU|) . and Eqs. (|6TT|) . ([BIZ]) , we have 



(*jv + ce) 



3.13) 



We are now ready to derive the differential equations 
for \ u so that (|6.2p are satisfied. We can cancel the 
common exp(i J2k ^k^kjR prefactor in (|6.10p and (|6.1ip 
(in the coordinate domain we are using the operator R 
is never zero on any spinor), multiply both equations by 

(— iy (— bi v \xj\ 1 3 and sum over j to obtain 



IXk 



/7-^-T/ \N-H>( X' $> v e \ 

{dxZ + 4t + X-^ + 2x-y 



2s?, 7 



'Xk • 



(6.14) 



Here we have used the 'completeness relations' 



E 



(0 



\ JV— 1-j ™ 



(.;) 



(6.15) 



together with the following identities: 

iV 

E* 

3=0 
N 

E B(J ' _1) (-w 

Using further the formula 



7 



we can rewrite Eq. (|6.14p as 



(6.16) 
(6.17) 

(6.18) 



dx v AX U X v ' s "' ' 2x v . 

?)X«=0. (6.19) 



2x1 



We finally note that the operators act only on the \v 
spinor, leaving invariant all the other spinors in the ten- 
sor product. So we are left with the following ordinary 
differential equation for each spinor \v'- 



f d X'* v ■ \ 

\d^^AX, + X^ + 2x-J 



\7V- 



(6.20) 



x, = o. 



This final equation is of the same form as for the un- 
charged case [lU; the only difference is the appereance 
of an extra term proportional to the coupling constant e 
in (HH) and fl6J]|. 

131 ] one can, by using the unitary trans- 



As discussed in 



formation which makes the 7 matrices coordinate depen- 
dent 



7 a = R-^ a R, 



3.21) 



achieve the standard separability (without prefactor R) 
in the form 



ij) = cxp(i J2k ^kipk 



3.22) 



where Xv satisfy the equation (|6.20p . The separated so- 
lution is an cigenfuction of operators Mj = R~ x MjR and 
Kk, i.e., a solution of 



(6.23) 
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To conclude, we have demonstrated complete separa- 
tion of variables for a charged Dirac equation in a weakly 
charged black hole spacetimes in all dimensions. Sep- 
arated equations (16.201) can be in principle decoupled, 
to yield a second order ordinary differential equation 
for each component of a 2-dimensional spinor Xv- It 
is easy to see that if instead of weakly charged black 
holes we considered weakly magnetized black holes the 
demonstrated procedure leading to the separated equa- 
tion (|6.20[) would break down at several places and hence 
the separability of the charged Dirac equation in this case 
is not possible, at least in the 'canonical setup'. We ex- 
pect this to be true in any reference frame. 



VII. SUMMARY 

Kerr-NUT-(A)dS rotating black holes in arbitrary di- 
mension are distinguished by the remarkable fact that 
several equations of physics in such a background can be 
solved by a separation of variables. This includes, among 
others, the Hamilton- Jacobi, Klein-Gordon and Dirac 
equations. Ultimately, the reason for this is the presence 
of a PKY tensor. It is possible to consider perturbations 
of these metrics by turning on an electromagnetic gauge 
potential along the directions of the Killing vectors. As it 
was shown in [l[ , among these only the perturbation that 
corresponds to a weakly charged black hole is such that 
it is compatible with separability of the Hamilton-Jacobi 
and Klein-Gordon equations. 

In this work we have studied the perturbed metrics 
from the point of view of the Dirac equation, analysing 
the presence of symmetry operators of the Dirac equa- 
tion and their relation to separability. A discussion of 
first order symmetry operators for the Dirac equation in 
the presence of an arbitrary flux can be found in (lo| . In 
particular, it is shown there how in order to build sym- 
metry operators for the Dirac equation in the presence 
of U{\) flux it is necessary to both i) have Killing- Yano 
tensors and ii) satisfy additional conditions of absence 
of anomalies. For all the perturbed Kerr-NUT metrics 
there is a candidate complete set of first order symme- 
try operators. However, we have proved that only for 
the weakly charged solution there are no anomalies. For 
this solution we have shown that all such operators mu- 
tually commute. We have also explicitly expressed the 
symmetry operators in a chosen frame, adapted to the 
hidden symmetry present in the spacetime, and demon- 
strated that in such a frame the Dirac equation admits a 
solution by separating variables. These results generalise 
those found in [l| and [H M, H • 
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Appendix A: Important identity 



In this Appendix we prove the identity p.22[) . We 
start with some auxiliary definitions. Let hr^) be the 
k-th "matrix" power of the PKY tensor h, 



(Al) 



For odd k, it is antisymmetric 2-form. (Note, is 
different from h^\) Since it is both Lie and parallel 
conserved along £, we have 

A h i2j+1) = . (A2) 

We also introduce 1-form with one vector index 

u,fa=X*_Lhv), (A3) 

which satisfies 

(A4) 



In the following we will need also two standard identi- 
ties for the contracted wedge product involving 1-form ui, 
p-form a and g-form /3: 

(lo Aa)A(3 = {-l) m uj A (a A (3) +ma A^wLi/3) , 

aA(w A/3) = (-lfuA (aA/3) +m(w'ja) A /3 . 

1 (A5) 



The identity (|3.19[) is based on the relation 

h (k) A h (i) = o (A6) 

2m+l V ' 

see the proof in [l6| . 

Now we can proceed in our proof of (|3.22p . First, we 
show the identity 

(V [o &]) wfo) A h^) A fcW) = . (A7) 

Expanding the mixed product with the help of (|A5[) . us- 
ing (|A4[) and the identity (|A6|) . we obtain 



(^L, h<») A fcW = u\ 2j+1) A (h^Ah^) 
+ (2m)(2m-l)( W( V 2) *j h^) a A_ a . 

Iterating this relation, we find 
(o^-ifcW) Ah® 



(A8) 



(2m)! 



(A9) 



S (2m-2i)\ 



-y (2]+2l+1) A (h^\ A_^) 
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Substituting back to (|A7|) wc obtain a sum of terms, each 
of which contains the expression^ 



(V[a&]) ^{2j) A UJ (2j+2 l +l) - h (2j) ■ (dt A fe(2i+l)) ' h(2j) , 

(A10) 

which vanishes thanks to relation (|A2[) . We thus proved 
the relation dATI). 



Next, we prove that 



(V [o &]) (e Q A 2 A +i (e 6 A ftW) = . (All) 



Employing again the relations (|A5|) for mixed products, 
the identity (|A6[) . the just proven relations (|A7|) . the an- 
tisymmetry of V[ a ^i,], and relation (|A4[) . we find 



(V [o fq)(e-AfcW)^ + > 6 A/i«) (A12) 
- (2m+l)(2m)(V [o &]) (wfo A M fc >) jvjwfo A . 



Defining e a = w" ), this expression suggests the recur- 
rence relation between ui^j) & n d ^(2j+2)i which can be 
proved by the same procedure. Iterating it (m + 1)- 
times, in the last step the expression annihilates com- 
pletely. 



Finally, the identity (|A1 1|) implies the desired relation 
p.22[) . Indeed, using the standard relations for Clifford 
anticommutator (see, e.g., (A19) of (la]), we get 



(V [o &]) [e^Ah^\e b ^Ah^]_ 
2(-l) 



(A13) 



£ (Lw ( v ^) ( e ° A - A > 6 A * w ) = 



Appendix B: Direct proof of mutual commutativity 
of operators Mj 



In this section we prove mutual commutativity of op- 
erators Mj, using their explicit representation (|5.16[) ; the 
presentation closely follows Sec. VI. A in [T3|. Let us start 
from the expression for Mj (|5 . 16[) and apply the identity 
(I6.18[) . together with 



E 



(Bl) 



to obtain 



.u. /'EvW'f 



ftc,, 4X M 2 



X 



f k 



B 



1 



et <M> 7" 



2 g 



(B3) 



1 



2e 



In order to prove commutativity of these operators we 
introduce new 'auxiliary' operators 

Mj = R _1 MjR , (B4) 

with R given by (|6.4[) . Then, obviously, if 

[M 3 -, M fc ] = R-^M^M^R = , (B5) 

then the same is true for the operators Mj . We calculate 

MjR = i^^/Q^B^RR^^R x (B6) 



<9:e u 4AT U 2x u X 



(7) 



1 



• f'f 



2 g 



-2e 



i?. 



where wc have used (|6.8j) . The i? factor in the e term on 
the right hand side can be brought to the front whereas 
for the product R~ 1 j fl R we use (|6.7[) . So we get 



B 



(./) 



Mj = i^^-M„ 



(B7) 



where the operators 



X„ 



dx^ 4 X M 2cc p 



et 



X 



* k 



— (- 

i /2 g 



</'> 



2e 



(B8) 



act only on spinor Xn an d hence [M^,M V ] = 0. Using 
(|6.15[) we can invert the relation (|B7[) . 

JV-l 

M M = ^ (-i)^-^^)^- 1 -^,- . (B9) 

Following now the procedure adopted in (29|, and using 
that [Mfj,, (~i.^x 1J ) N ~ 1 ~ : >] = 0, we establish that 

N-l 

x; (- J ) j+fe (-^ ) ^) jv - i --'"(- t< , )a; ,) Ar - i - fc x 

j,k=0 

x[Mj,M k }=0, (BIO) 

from which Eq. (|B5|) follows. 
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